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Within the framework of R-parity violating minimal supergravity model, at least three relatively 
large lepton-number violating A' type trilinear couplings at the GUT scale, not directly related 
to neutrino physics, can naturally generate via renormalization group (RG) evolution and/or CKM 
rotation the highly suppressed bilinear and trilinear parameters at the weak scale required to explain 
the neutrino oscillation data. The structure of the RG equations and the CKM matrix restrict the 
choices of the three input couplings to only eight possible combinations, each with its own distinctive 
experimental signature. The relatively large input couplings may lead to spectacular low energy 
signatures like rare weak decays of the r lepton and K mesons, direct lepton number violating 
decays of several sparticles, and unconventional decay modes (and reduced lifetime) of the lightest 
neutralino, assumed to be the lightest supersymmetric particle (LSP), all with sizable branching 
ratios. Several low background signals at the Tevatron and LHC have been suggested and their 
sizes are estimated to be at the observable level. From the particle content of the signal and the 
relative rate of different final states the input couplings at the GUT scale, i.e., the origin of neutrino 
masses and mixing angles, can be identified. 

PACS numbers: 12.60.Jv, 14.60.Pq, 13.85.Rm 



I. INTRODUCTION 

The discovery of neutrino oscillations [1] has estab- 
lished beyond doubt that there are at least two massive 
neutrinos. However, their masses are much smaller than 
the masses of the other known fermions. Discovering the 
origin of these small masses is perhaps the most challeng- 
ing task for current high energy physics. 

The see-saw mechanism which can be naturally acco- 
modated in any grand unified theory (GUT) is one of 
the most popular explanations of small v masses [2] . Un- 
fortunately the simplest version of this theory, a GUT 
with a grand desert, has very few predictions for low 
energy physics apart from v masses and mixing angles. 
In particular the low energy spectrum of such a theory 
is practically identical with that of the Standard Model 
(SM). 

Supersymmetry (SUSY) is the most elegant extension 
of the SM which solves the naturalness problem that is 
inevitable in any non-supersymmetric GUT [3]. There 
is a version of the minimal supersymmetric extension of 
the SM (MSSM) which conserves lepton number L and 
baryon number B due to an additional discrete symmc- 
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try called R-parity and defined in such a way that all 
particles have R= 1 and all sparticles have R= —1. This 
theory is referred to as the R-parity conserving (RPC) 
SUSY. However, there are other interesting variations of 
the MSSM with appropriate discrete symmetries which 
violate either baryon number or lepton number [4] but 
not both (so proton decay does not occur). These vari- 
ants are known as the R-parity violating (RPV) SUSY. 
In the most general version of this theory there are bilin- 
ear and trilinear RPV terms [4] in addition to the usual 
RPC terms. 

One of the most interesting features of the L- violating 
version of the RPV SUSY is that it can naturally ex- 
plain the experimentally measured v masses and mixing 
angles [5]. More important, as demanded by the natu- 
ralness argument, this model is entirely governed by TeV 
scale physics. Thus the particle spectrum consists of su- 
perpartners, i.e., the sparticles, of the particles of the 
SM, having masses in the TeV scale. The production 
and detection of these sparticles at the ongoing Tevatron 
run II and the upcoming Large Hadron Collider (LHC) 
or the International Linear Collider (ILC) accelerator ex- 
periments can directly test the RPV models of v mass. 

In the RPC SUSY the lightest supersymmetric particle 
(LSP) is stable. In contrast RPV SUSY allows the LSP to 
decay, through the RPV couplings, to lepton number vio- 
lating channels. The multiplicity of particles in any event 
involving sparticle production is, therefore, much larger 
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on the average. The other characteristic signatures [4] of 
RPV SUSY are single production of sparticles, direct de- 
cays of sparticles via lepton number violating channels, 
and in particular the decay of the lighter top squark [4, 6- 
9]. 

In addition to the above direct tests, indirect signa- 
tures of RPV or RPC SUSY can be obtained through 
rare weak processes such as highly suppressed K, B, D, 
or r decays or decays forbidden in the SM [4, 10, 11]. 
The standard signal is an abnormal enhancement of 
the branching ratio (BR) and/or significant change of 
CP asymmetries compared to the SM expectations [12]. 
Moreover the observed K° — K° and B° — B° mixings 
put strong constraints on RPV parameters [13-15] as we 
shall see later. It should be emphasized that RPC SUSY 
can contribute to these processes only at the loop level, 
while RPV SUSY can contribute at the tree level. The 
latter contributions can, therefore, be potentially large 
and predict significant deviation of various observables 
from the corresponding predictions in the SM. 

No evidence of RPV or RPC SUSY has so far been 
found either through direct or indirect methods. This 
leads to bounds on the parameter space of RPV SUSY 
[4, 16]. Particularly interesting are the bounds on the 
trilinear RPV couplings which control the size of various 
direct or indirect signatures. Some of these bounds, de- 
rived before the advent of the highly constrained neutrino 
oscillation data [1] or of precise cosmological observations 
[17], were remarkably weak and the corresponding RPV 
couplings could be as large as allowed by the perturba- 
tive nature of the theory! The magnitudes of direct or 
indirect signatures of RPV SUSY based on these bounds 
were, therefore, overestimated. 

Many authors have revisited the bounds in the light 
of the neutrino data [18-20]. A major problem of the 
most general RPV model is that the number of free pa- 
rameters are even larger than the almost unmanagably 
big parameter space of the RPC MSSM. In order to make 
the analysis tractable additional simplifying assumptions 
restricting the parameter space of the RPV sector were 
employed. 

One approach is to consider a few selected benchmark 
scenarios of v oscillation. Each scenario consists of a 
minimal set of RPV bilinear and trilinear couplings at the 
weak scale, the number of parameters being just enough 
to cope with the oscillation data [18]. Upper bounds on 
the parameters belonging to each set were then obtained 
by a global analysis of the v oscillation data. 

In this paper we shall focus on the scenario with three 
bilinear m and three L-violating trilinear A - 33 couplings 
(the leptonic index i will run from 1 to 3 in the entire 
paper). While the bilinears generate a tree- level mass 
matrix for the neutrinos, the trilinears generate such a 
matrix at the one loop level. The upper bounds on the 
A' couplings turn out to be rather strong (~ 10~ 4 ) [18]. 
As a result the contributions of these couplings to most 
low energy processes except LSP decay are negligible. 
One notable exception is the direct RPV decay of the 



lighter top squark [8, 9] if it happens to be the next to 
lightest supersymmetric particle (NLSP), a theoretically 
well-motivated scenario. This is so because the compet- 
ing RPC decays of the top squark NLSP are also natu- 
rally suppressed (see Section V for the details). Of course 
other RPV phenomena at observable levels can be accom- 
modated in ad hoc models by arbitrarily adding other 
relatively large couplings not directly related to v phe- 
nomenology Such phenomena, however, are not corre- 
lated in any way with the v sector and, hence, can throw 
no light on the origin of v masses and mixing angles. 

Another attractive approach is to assume that the fun- 
damental theory at some high scale, say the GUT scale 
Mq, has a small number of parameters, of RPV or RPC 
type. A larger set of parameters required by low en- 
ergy phenomenology, including v oscillation, can be gen- 
erated at the weak scale via the renormalization group 
(RG) evolution of these parameters. One popular way 
is to consider the usual boundary conditions at Mq of 
the minimal supergravity (mSUGRA) model along with 
three bilinear RPV parameters only [21]. However, these 
models have very little predictive power beyond LSP de- 
cay which has been studied in detail [22]. 

Subsequently the full set of both bilinear and trilinear 
RPV couplings were also included in the analysis [23]. 
A recent work [20] have obtained strong bounds on A' 
couplings from the WMAP [17] constraints on the sum 
of v masses. They extended the conventional mSUGRA 
boundary conditions at Mq minimally by adding only 
one non-zero A' couping. 

It should be noted that the conventional GUTs do not 
include gravity. Thus any GUT may be regarded as an 
effective theory embedded in a more fundamental theory 
of gravity. It is quite conceivable that this theory will 
lead to non-renormalizable opoerators at Mq suppressed 
by a heavy mass scale [24]. When the GUT symme- 
try is broken down and some heavy Higgs fields develop 
vacuum expectation values (VEV), bilinear and/or tri- 
linear RPV interactions are generated at Mq- The mag- 
nitudes of these parameters depend, among other things, 
on these VEVs. Due to our rather limited knowledge of 
the high scale physics at the moment the magnitudes of 
the induced RPV parameters cannot be computed from 
first principles. However, it is quite possible that only a 
small subsets of these parameters have numerically signif- 
icant magnitudes. At the moment a meaningful question 
would be to ask what is the smallest set of RPV param- 
eters which can account for low energy phenomenology 
including neutrino physics. The magnitudes of these pa- 
rameters can be restricted by low energy data. 

In the model of ref. [20] the single A' coupling at Mq 
induce RPV bilinears at the weak scale via RG evolution. 
This in turn induces one tree level v mass at the weak 
scale incompatible with the WMAP bound unless the 
input A' coupling is very tightly bounded from above. It 
turns out that almost all the 27 A' couplings have tiny 
upper bounds [20] so that the predictions for the indirect 
signatures of RPV and/or the direct RPV decays of the 
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sparticles, except for the LSP, are not at observable levels 
attainable in the foreseeable future. 

However, a single A'-type coupling cannot generate the 
full 3x3 neutrino mass matrix at the weak scale with 
contributions from both tree and one loop amplitudes, 
which could in principle be of the same order of magni- 
tude. Thus the stringent bounds of [20], though highly 
suggestive, should not be regarded as a general conclu- 
sion following from the RPV models of m„. 

In this paper we propose a novel mechanism with a 
minimal set of relatively large A' couplings, not directly 
related to v masses, as inputs at Mq. Thanks to flavor vi- 
olating effects like non-diagonal Yukawa coupling matri- 
ces, this minimal set at Mq can induce, via RG evolution, 
several new non-zero A' couplings in the weak or flavor 
basis at the weak scale with naturally suppressed mag- 
nitudes. The couplings in the physical or mass basis are 
obtained by applying appropriate Cabibbo-Kobayashi- 
Maskawa (CKM) rotations. In addition several RPV bi- 
linear terms are generated, even though they are zero at 
Mq according to the chosen boundary conditions. While 
some of the naturally suppressed induced RPV param- 
eters can take care of the v oscillation data, some rela- 
tively large A' couplings corresponding to the inputs at 
Mq exist at the weak scale with spectacular testable con- 
sequences for low energy physics. 

It should, however, be emphasized that the choices of 
the minimal input sets are indeed very much restricted. 
This will be explained in Section III by analyzing the 
structures of the RG equations and the CKM matrix and 
by taking into account the experimental constraints cur- 
rently available. 

As an illustrative example we focus on the benchmark 
scenario consisting of three bilinear parameters Kj, and 
three trilinear couplings A^ 33 [18] at the weak scale. We 
find that the set of input couplings for inducing this 
benchmark scenario should be of the form A^ 3 , i = 1, 2, 3 
and j = 1,2. Any three of these six input couplings bear- 
ing different lepton indices provide the minimal set we 
wish to find out. 

For the purpose of illustration we have focussed on a 
particular benchmark scenario of v oscillation. One can 
in principle start with other appropriately chosen mini- 
mal set of input parameters {e.g., three A type couplings) 
at Mq having relatively large magnitudes and induce the 
the desired benchmark scenario (i.e., Ki and Ai33, [18]) at 
the weak scale if CKM like flavor violation in the lepton 
sector is taken into account. Of course v phenomenol- 
ogy alone cannot distinguish between various benchmark 
scenarios, but the phenomenology of the input couplings, 
which we shall discuss next, automatically points to the 
underlying model of v oscillations. 

The magnitudes of the input couplings can in principle 
be determined by the data on v masses and mixing an- 
gles [1]. Unfortunately the detailed predictions for the v 
sector depend also on the parameters of the RPC sector 
(see below for the details) . Practically no information on 
these parameters are available at the moment. 



We hope that sufficient information will soon be avail- 
able on the RPC sector from the measurement of sparticle 
masses and BRs at Tevatron run II and LHC. Moreover, 
v data will become more precise in course of time. Only 
then a complete phenomenological fit of the model pa- 
rameters to the oscillation data can be possible. For the 
time being, in order to estimate the size of the signals 
predicted by our model, mainly for the purpose of illus- 
tration, we have restricted the input couplings in a rough 
way. We require that the magnitudes of the input cou- 
plings at Mq be such that the induced parameters at the 
weak scale related to the neutrino sector satisfy the up- 
per bounds derived in [18] from oscillation data. More 
discussions on this point can be found in Sections II and 
V. 

The main thrust of this paper is to identify the mini- 
mal sets of input couplings at Mq which are still allowed 
to be relatively large in spite of the severe constraints 
from v oscillation data. We then examine the interesting 
rare weak decays and collider signatures, including direct 
RPV decays of sparticles other than the LSP or the top 
squark, which are allowed to be at the observable level 
in our model. 

In general the elements of the Yukawa coupling matri- 
ces in the up and the down quark sector, denoted by Y u 
and Yd respectively, cannot be determined from the quark 
masses and the measured magnitudes of the elements of 
the CKM elements. On the other hand, these matrices 
appear explicitly in the RG equations. Quite often ad- 
ditional simplifying assumptions on the structures of the 
quark mass matrices are introduced so that the elements 
of Y u and Yd are calculable from the observables. 

In the first part of our analysis, summarized in the 
above paragraphs, we have assumed that the quark mass 
matrices have a particular texture such that the CKM 
matrix is identical to the mixing matrix in the up quark 
sector [20, 25]. If on the other hand a scenario with 
mixing restricted to the down quark sector is considered, 
the magnitudes of the input couplings at Mq are very 
severely restricted by the bounds of [18] and can hardly 
lead to any observable low energy phenomenology apart 
from the LSP decay. This feature of the second scenario 
has also been noted in [20], though their conclusion was 
based on a model with a single neutrino mass. 

We then take a more phenomenological approach and 
assume that some RPV bilinears as well as A' couplings 
in the weak basis, not directly related to v physics, are 
generated at low energies with relatively large magni- 
tudes by some high scale physics which is not necessarily 
of supergravity type. Starting from a limited number of 
input coupling sets at the weak scale, which are different 
in general from the ones obtained in the first scenario, 
the parameters required to explain the oscillation data 
can in principle be generated through CKM rotations. 
However, such rotations will generate other parameters 
as well. Using the currently available constraints on the 
elements of the CKM matrix [26] and those on the prod- 
uct RPV couplings arising from K° — K° and B° — B° 
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mixing, we find that the input couplings are so severely 
restricted that neither the rare decays nor the direct lep- 
ton number violating decays of sparticles (other than the 
decay of the top squark or the LSP) triggered by them 
will be observable. However, both the top squark and the 
LSP decays will have distinctive features characteristic of 
the second scenario which will be discussed in detail. 

The plan of the paper is as follows. In Section II we 
establish the notation, and discuss briefly the method- 
ology of solving the relevant RG equations and determi- 
nation of the physical RPV couplings from the ones in 
the weak basis generated by RG evolution in Scenario I 
(mixing restricted to the up quark sector). In Section III 
we identify the relatively large input trilinear RPV cou- 
plings at Mq in Scenario I which can induce at the weak 
scale the naturally suppressed parameters for neutrino 
oscillations and other interesting couplings via RG evo- 
lution and CKM rotation. In Section IV we discuss the 
testable predictions of our model for rare decays. In Sec- 
tion V the RPV decays of the top squark, the LSP, and 
other sparticles, and their collider signatures have been 
analyzed. In Section VI the phenomenological model at 
the weak scale in Scenario II (mixing restricted to the 
down quark sector) and its observable consequences are 
examined. Our conclusions are summarized in the last 
Section. 



II. RG EVOLUTION AND CKM ROTATION OF 
THE RPV COUPLINGS 

With an eye on the neutrino phenomenology, we con- 
fine ourselves to L-violating couplings only. There are 
three such terms in the superpotential: 

W 9 = hjkULjEl + \' ijk LiQjDl + (1) 

where L, E, Q, D and H 2 stand for doublet lepton, singlet 
charged lepton, doublet quark, singlet down-type quark, 
and the second Higgs doublet (that gives mass to the 
top quark) superfields, and i, j, k = 1, 2, 3 are generation 
indices. The first term is antisymmetric in i and j. 

Regarding the RG evolution of the RPV couplings, we 
follow [20] and refer the reader to this paper for the de- 
tails. The anomalous dimensions of the standard Yukawa 
couplings as well as of the RPV couplings, and the (3- 
functions, including that of the bilinear terms fi and Ki, 
are given in the appendix of [20] , and we will not repeat 
them here. 



For the benchmark scenarios at the weak scale we are 
interested in, it is sufficient to consider only the A^ fc type 
couplings at Mq- It is obvious from the RG equations 
that A^jfcS are never generated from A- jfc s through RG 
equations and vice versa. However, even if the bilinear 
couplings Ki are zero at the GUT scale, they will be gen- 
erated at the weak scale mainly through a term in the 
RG equation that depends on the trilinear \[- k couplings. 

The first part of the procedure is to find the RPV pa- 
rameters at the weak scale. This is done through a mod- 
ified version of the code ISAJET v7.69 [27] to which we 
have added the RG equations for the RPV parameters. 
The steps are outlined below in brief. 

• The gauge and Yukawa couplings are specified at 
the weak scale. We assume that the leptonic 
Yukawa matrix is diagonal. As discussed in the in- 
troduction some simplifying assumptions must be 
made about the structures of the Yukawa coupling 
matrices Y u and Yd in the quark sector so that their 
elements, which appear in the RG equations, can 
be directly determined from the quark masses and 
mixing angles. The quark Yukawa matrices at the 
weak scale are assumed to be real and symmetric 
so that the rotation matrices in left and right sec- 
tors are the same. In the most general case The 
Cabibbo-Kobayashi-Maskawa (CKM) matrix V is 
given by V = U' l Dl where Ul and Dl are the rota- 
tion matrices for left-handed quark fields. Follow- 
ing earlier works [20, 25] we consider two scenarios: 
Scenario I: U L = U R = V\ D L = D R = 1; 
Scenario II: D L = D R = V, U L = U R = 1. 

Subject to the above assumptions the entries of the 
Yukawa matrices can be computed from 

(7) m d (M z ) = [m d ] diag (M z ) , 

m„(M z ) = V* KM • [m u ] diag (M z ) • V£ KM . 
(II) m d (M z ) = V* CKM ■ [m d ] diag (M z ) • V% KM , [ > 

m u (M z ) = [m„] diag (M z ) . 

For simplicity, we neglect u and d quark masses. 
One can also neglect the second generation quark 
masses, and all lepton masses except that of r. The 
magnitude of the entries of the CKM matrix are 
taken from Particle Data Group 2004 [26], based 
on the unitarity of the matrix: 



'0.9739- 0.9751 0.221 - 0.227 0.0029 - 0.0045\ 
Vckm=\ 0.221 - 0.227 0.9730 - 0.9744 0.039 - 0.044 . (3) 
0.0048 - 0.014 0.037- 0.043 0.999 - 0.9992 J 

I 



For numerical calculations we consider the central 



value of each element, except V U b and V t d- 



■5 



• At the GUT scale the elements of Y u and Yd can 
be computed from the low energy inputs (eqs. (2) 
and (3)) via the RG equations following an itera- 
tive procedure [20] . The elements so computed will 
govern the RG evolutions of the RPV parameters. 
The numerical values of Y u {Mg) and Yd(Mc) in 
Scenario I will be presented in Section III. 

Scenario II will be discussed in Section VI. 

• The mSUGRA parameter space is specified by mo, 
the common scalar mass, m 1 / 2 , the common gaug- 
ino mass, tan/3, the common trilinear term Aq, and 
sgn(p). We run all the gauge and Yukawa cou- 
plings upwards until g\ and g 2 meet; that is taken 
to be the unification point. The strong coupling 
g 3 is unified by hand at that point. It is checked 
whether gz{M z ) in turn goes outside the experi- 
mental range; it does not. The mSUGRA parame- 
ters used for the computation of various signals will 
be presented in Section V. 

• At the GUT scale, all mSUGRA parameters are 
specified. We also specify some A' type RPV cou- 
plings to be nonzero at Mq- As we have mentioned 
earlier and shall further clarify in the next section, 
we need to specify at least three nonzero A's each 
with a different lepton index at the GUT scale and 
in the weak eigenbasis of the quark fields in order 
to successfully reproduce the v phenomenology at 
the weak scale. 

• We run down to the weak scale, taking all thresh- 
old effects into account. It is ensured that the 
electroweak symmetry is radiatively broken deter- 
mining thereby the magnitude of /x. We also per- 
form some sample checks on the scalar directions of 
the superpotenial at various energy scales so that 
no dangerous directions with unbounded or charge- 
color breaking minima are encountered. Of course, 
this is not a serious objection if we live in a false 
vacuum, but the solution is not aesthetically pleas- 
ing. 

• The above running of the input RPV couplings in- 
duce other trilinear and bilinear couplings at the 
weak scale including m and A- 33 required for the 
neutrino sector. 

• One notes the crucial role the quark mixing matri- 
ces play: without them, no A^ fc which is zero at the 
GUT scale can acquire a nonzero value at the weak 
scale. This can easily be checked from the RG equa- 
tions as given in [20]. Thus, with the CKM matrix, 
our minimal set of 3 A' couplings at the GUT scale 
can generate all 27 couplings to be nonzero at the 
weak scale. These couplings are in the weak eigen- 
basis and should be rotated by the proper quark 
mixing matrix to get the physical couplings in the 
mass eigenbasis. How this is done is shown in the 
next section. 



• The RG evolution of the input A' couplings are by 
and large independent of the RPC sector of the 
mSUGRA parameters apart from tan (3 which re- 
lates the quark masses and the Yukawa couplings 
(eq. 2). The magnitudes of the KiS at the weak 
scale, however, also depend on the higgsino mass 
parameter /x of the RPC sector. Usually the mag- 
nitude of this parameter is fixed by the radia- 
tive electroweak symmetry breaking condition [28] . 
This magnitude depends, among other things, on 
the assumption of a common soft breaking scalar 
mass mo at Mq. There are ample reasons to be- 
lieve that even within the supergravity framework 
the soft breaking masses for the Higgs sector and 
the sfermion sector may be significantly different 
[29, 30] at Mq. In such nonuniversal models the 
magnitude of \x may be quite different from that in 
the mSUGRA model. This in turn introduces siz- 
able uncertainties in the magnitudes of the bilinear 
RPV parameters. Moreover, both the tree level 
and one loop neutrino mass matrices depend on 
the mSUGRA input parameters which are poorly 
known at the moment. 

• We have studied the variation of X' i j k and K{ with 
tan/3. If tan/3 is changed from 5 to 20, the cou- 
plings which are important in the study here, viz., 
X' il3 , X' i2 3 and A^ 33 , change by less than 10%. Hence 
our estimates of the sizes of various collider signals 
(Section V) remain practically unaffected. How- 
ever, (Yd)ij at Mq changes significantly, and as a 
result there is an appreciable change in Kj. Thus if 
one attempts to fit the v oscillation data precisely, 
tan (3 dependence must be taken into account. 

• We do not attempt a detailed RPV parameter fit- 
ting using the v oscillation data for reasons dis- 
cussed in the last paragraph. We, however, ensure 
that the input RPV parameters have the right ball- 
park values by satisfying the following criteria: 

(i) The bilinear couplings KiS at the weak scale, 
related to the elements of the tree level neutrino 
mass matrix, should not violate the upper limits 
(~ 10- 4 GeV) [18]; 

(ii) The magnitude of the trilinear couplings A^ 33 , 
which are responsible for the one loop neutrino 
mass matrices, should also be bounded from above 
(< 1.5 x 10~ 4 ), so that the neutrino mass squared 
splittings and the mixing angles are reproduced. 

Strictly speaking the bounds of [18] were derived for 
a common SUSY scale M SUSY = 100 GeV for all 
masses and mass parameters and should be mod- 
ified when an mSUGRA mass spectrum is consid- 
ered. But since we are mainly interested in right 
ballpark values only, we ignore possible changes in 
the bounds, which are not expected to be very dras- 
tic. Moreover, the bounds of [18] were derived in a 
basis in which the sneutrino VEVs vanish, whereas 
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after the RG evolution from the GUT scale to the 
weak scale we arrive at non-zero values of the above 
VEVs in general. In principle a rotation in the 
ifd-slepton space should be applied to make these 
VEVs zero. However, such rotations are expected 
to change the couplings by factors which are order 
one. We neglect such effects in our order of magni- 
tude estimates. 

(iii) In the mass eigenbasis, all the existing con- 
straints on the individual couplings and their prod- 
ucts [4, 16] from non-neutrino physics should be 
satisfied; 

(iv) Clearly the v oscillation data cannot be ex- 
plained if all the induced parameters and A- 33 ) 
have magnitudes much smaller than the above up- 
per bounds. In order to estimate roughly the 
magnitudes of the input couplings, we have as- 
sumed that the magnitudes of the above param- 
eters should be of the same order of magnitude 
of the upper limits mentioned in (i) and (ii). In 
a forthcoming paper [31], based on numerical cal- 
culation of the eigenvalues and eigenvectors of the 
neutrino mass matrix (including both tree level and 
one loop contributions) and comparison of the re- 
sults with the measured v mass squared differences 
and mixing angles, it will be shown that the above 
assumption is fairly realistic (also see Section V). 

Let us now discuss how far the v data can be useful 
in constraining the RPV scenarios under consideration. 
First take Scenario I. In the mass eigenbasis, eq. (1) reads 

W 9 = -M ,/, A , / ), /); + % mk EiU m D° p , (4) 

with 

\/ \/ T"' \/ ta* /r\ 

ijk ijk' "iink ^ijk* jm' \P ) 

In eq. (4) Ni and Ei are respectively the neutral and 
charged components of the SU(2) doublet lepton super- 
field Li. Thus, in the ^-sector, new couplings must be 
generated from the inputs at Mq through RG evolution 
alone; CKM rotations do not play any role. However, 
both RG evolution and CKM rotations may combine to 
generate at the weak scale interesting couplings in the 
mass basis involving charged lepton fields (see Sections 
IV and V for the consequences). Also see Section VI for 
Scenario II. 



III. CHOICE OF THE INPUT RPV 
PARAMETERS AT M G 

Our goal is to find the minimal set of trilincar RPV 
couplings at Mq capable of reproducing at the weak scale 



the benchmark scenario [18] for v oscillation under con- 
sideration. 



Substituting for the anomalous dimension matrices in 
the RGE's of the trilincar couplings, eq. (A23) of [20], it is 
obvious that in order to generate a particular A^ fe at the 
weak scale we need a nonzero X' imn (with the same lepton 
index) at Mq. Thus at least three A'-type couplings each 
bearing a different lepton index are required as inputs at 
Mq. The same conclusion follows if we examine the RG 
equations of the KiS, and require these parameters to be 
non-zero at the weak scale for all values of i. 



Of course, one option is to consider three couplings 
A^ 33 as the inputs at the GUT scale. However, the mag- 
nitudes of these inputs will be so small due to neutrino 
constraints that no new couplings having magnitudes ap- 
propriate for observable signals at the weak scale will be 
induced. Thus the only collider signals will be LSP de- 
cays [32] and top squark decays [8, 9] triggered by the 
three A^ 33 -type couplings themselves. These have, how- 
ever, been studied in detail and we shall not consider this 
option further in this paper. 

Thus the desired minimal set of input parameters at 
Mq may contain any three of the 24 A^ fe -type couplings 
(with j = k = 3 excluded) each with a different lepton 
index. First of all we reject A' in because from neutrino- 
less double (3 decay we get a very strong upper bound on 
this (< 0(1(T 5 )). 

The Yukawa matrices Y u and Yd at Mq are needed to 
study the RG equations of the RPV parameters. Follow- 
ing [20] the elements of the matrix Y u (Mq) are deter- 
mined from the quark masses and mixing angles. Their 
magnitudes in Scenario I (eq. (2)) are found to be 



/l.67xl0~ 4 7.4 x 10~ 4 8.4 x 10~ 5 \ 
Y U (M G ) = 7.36 x 10~ 4 3.5 x 10~ 3 2.3 x 10~ 2 , 
\-9.8 x 10~ 6 2.3 x 10~ 3 5.7 x lO" 1 / 

(6) 

where we have used m u = rrid ~ 0.0, m s = 0.199 GeV, 
m c = 1.35 GeV, TO fc = 4.83 GeV and m t = 175 GcV. The 
off-diagonal elements of Yd at Mq are generated through 
mixing in the up sector as seen from the equations given 
below. Starting from a diagonal Yd at the weak scale we 
have 
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i67r 2 |(r d ) 22 = (Y d ) 22 

l67r 2 j t (Y d ) 33 = (Y d ) 33 



• :!//;: • ) • - -.>ir, ): i, - o;.);,, 

^ + 3s 2 2 + ypl) + |(r„) :! „| 2 t--.V/M, MH>;/)r;, + (i; )l- ( 



(7) 



For the off-diagonal elements of Y d the equations are 

These three equations generate the elements of Y d (Ma). 
The only numerically significant elements are (Id) 13 = 
9.5 x fCT 6 , (Y d ) 22 = 3.5 x f(T 3 , (Y d ) 23 = 2.5 x f(T 4 , 
{Y d ) 32 = 1.73 x 1(T 5 , and (Y d ) 33 = 5.64 x 1CT 2 ; all other 
elements of Y d are negligible. 

The procedure to identify the phenomenologically in- 
teresting minimal input sets is now outlined below. 

• As discussed in the last section, we require A- 33 
and Ki to be of the same order of magnitude of 
their upper bounds, 1.5 x 1(T 4 and 1.2 x 1CT 4 GeV 
respectively at the weak scale. 

• If we consider the X' 211 and X' 311 couplings as inputs 
at Mq then from the RG equation 

16^A; 33 = -3 x (FcOssTOii^u (9) 

it follows that the evolution of A- 33 is controlled by 
the (Y d )n which is too small to give A- 33 with the 
desired order of magnitude. That is why we ignore 
A 211 and A 311 as possible inputs. 

• If we numerically integrate the RG equations for 
A^ 33 and n t using any of A^ 21 , A- 31 or \' il2 as inputs, 
we find that the last condition can be satified only 
if the input couplings are very large and grossly 
violate the existing upper bounds [4, 16] from non- 
neutrino phenomenology. This can be understood 
by the following qualitative argument. 

Suppose we integrate the RGE's of ^ using the 
crude approximation that all couplings other than 
Ki are constants having their respective values at 
Mq at all scales. We obtain 

Ki w = -l.2^f{{Y d ) nm ) G {\' mm ) G . (10) 

where the superscripts W and G refer to the weak 
scale and the GUT scale respectively. The con- 
stant 1.25 contains the large logarithm involving 
M G = 2 x 10 16 GeV, M z = 91.1 GcV and other 
multiplicative constants appearing in the RG equa- 
tion [20]. Substituting the values of (^) G and 
{(Xd) nm ) G and using m as restricted above, we es- 
timate the values of X' inm (with n = m = 3 ex- 
cluded) at Mq. We find that due to small values of 



the corresponding elements of Y d (Mc), the param- 
eters A- 31 , A- 12 and A^ 21 become too large. On the 
other hand \' i22 at Mq is required to be too small 
to be of any phenomenological interest, due to the 
relatively large value of ((Y d ) 22 ) G . 

• The remaining choices are A- ;3 (I = 1,2) and A^ 32 . 
Using cq. (10) we can estimate the values of the 
couplings A^ 3 and A^ 32 , which induces kY in the 
right ballpark. Integrating the RG equation for A^ 33 
using the same approximation as above, we find 

\' l33 W ~ 2.09A^ 32 G x (Y d ) 33 G x (F d ) 32 G . (11) 

When the estimated values of A^ 32 are substi- 
tuted, X' i33 W turns out to be O(10~ 7 ) which is too 
small to be of any interest in neutrino physics. On 
the other hand, if the elements (A^ 3 ) G are nonzero 
at M<3, a rough estimate of A- 33 at the weak scale 
is 

A^ ~ 0.42A^ 3 G x (Y u ) 3n G x (Y u ) ln G . (12) 

Compairing magnitudes of the elements of the ma- 
trices Y u and Y d at Mq, we conclude that the es- 
timated values of A^ 3 naturally reproduce values 
of A^ 33 correct to the order of magnitude estimate 
that is favored by v data. 

So, at the end of the day, we have only six couplings, 
A- 13 and A- 23 , out of which any three with different lepton 
indices constitute a minimal set of inputs at Mq. From 
the numerical solutions of the RG equations we find that 
|A- 13 | < 0.13 and |A^ 23 | < 0.26 x 10~ 2 can induce k, and 
•M33 w ith the correct order of magnitude at the weak 
scale. 



IV. RARE WEAK DECAYS 

In order to study the rare weak decays in Scenario I, 
we start with a set of three A'-type couplings in the weak 
basis as inputs at the GUT scale (Section III). For suc- 
cessful explanation of v phenomenology, these three cou- 
plings should have different lepton indices. After running 
down to the weak scale, as has been outlined in Section II, 
we get nonzero values for all A'-type couplings, but they 
are still in the weak basis. For the interaction involving 
v fields no further CKM rotation is required to get the 
interactions in the mass basis. For the charged lepton 
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interactions, on the other hand, we rotate the fields to 
the mass basis (eqs. (4) and (5)) using the CKM ele- 
ments given in eq. (3). The initial values of the three 
nonzero couplings at Mq are so chosen that at the weak 
scale all |A^ 33 < 1.5 x 10 -4 . We start with the largest 
possible values thus allowed, and check whether any in- 
dividual induced coupling X imk at the weak scale and in 
the mass basis is in conflict with the already existing ex- 
perimental upper bound[4, 16]. If that is the case, we 
scale down the input couplings accordingly. Otherwise 
we check whether any product of two such A couplings 
violates the corresponding bound, and if that is so, we 
scale down the input couplings further. This iterative 
process is continued till all couplings at the weak scale 
(and their products too) are consistent with the exper- 
imental limits. We next study a few specific choices of 
inputs at Mq. 

Case 1: AJ 13 

For this case, to start with, \X' il3 (Mc)\ could be at 
most 0.13 due to oscillation constraints alone. Their 
magnitudes increase about threefold at the weak scale. 
However, even the individual bounds on A' 113 and A' 213 
are tighter provided the squarks arc not too heavy. 
The first one is 0.02(m^/100) and the second one is 
0.06(mg fi /100) [4, 16] where is the generic sfermion 
mass. The bound on A' 113 comes from charge current 
universality (CCU) (as well as atomic parity violation 
(APV)) and that on A' 213 comes from the measured ra- 
tio R„ = r(7r -> ev)/r(ir -> \iv) [4, 16]. After taking 
this scaling down into account, the RG induced couplings 
| A^ 33 1 relevant for the v sector now have acceptable order 
of magnitude values. 

The charged lepton couplings |A- 13 | in the weak 

basis change to the A couplings in the mass ba- 
sis due to CKM rotation. The new couplings thus 
generated violate the upper bounds on the products 
1*313 AsssKiOO/m,) 2 < 2.0 x 10~ 3 coming from B° - B° 
mixing and |A 313 A 323 |(100/to ; ~) 2 < 2.7 x 10~ 3 from 
K° — K° mixing [15] unless the input values arc fur- 
ther scaled down (for the relevent formulae see eq.(21) 
of [15]). Only these constraints, and not the ones from v 
oscillation, control the ultimate upper limits of the input 
values of the A's at the GUT scale. Our final choice is 
A' 113 (100/m 5 J < 0.0064, A 213 (100/m^ R ) < 0.0194, and 
A 313 (100/m r ) < 0.0386. 

Such a set of couplings generate several novel channels 
for rare r decays, e.g., r — > pp, which is forbidden in the 
SM. The decay amplitudes can be found in cq. (11) of 
[15]. The BR scales with m^ 4 . In Table (I) we show the 
maximum theoretical expectations vis-a-vis the experi- 
mental upper bounds [26] for a number of lepton flavor 
violating channels, with inputs mg H = 300 GeV and rrif 
= 100 GeV. It is interesting to note that the theoretical 
limits turn out to be only one order of magnitude (or 
even less) smaller than the experimental data in many 
cases. 



J. I l_l l_> L 

coupling 


T Tnnpr bound 
at M w 


A 1 VJ V V .111 


BR 


Expt. 
limit 


^313^113 


2.4x 10 -3 


r — > en 

T ► CTj 
T — > Cp 


8.5xl(T B 
7 a v 1 n -8 

I .OA 1\J 

1.5xl0" 7 


3.7xl(T b 
8.2x 10 -6 
2.0xl0 -6 


^313-^213 


7.2xl(T 3 


t — ► pn 
t -> pr\ 
T * MP 


7.8x10"' 
7.2xl0 -7 
1.4xl0 -6 


4.0xl0 _b 
9.6xl(T 6 
6.3xl0 -6 



TABLE I: RPV mediated rare decays of the r lepton with 
possibly large BRs. Since A^ 13 is the input coupling set, there 
is hardly any difference between A' and A for these couplings. 



The trilinear couplings of ^s with two down-type 
quarks do not feel any effect of the CKM rotation, and the 
nonzero values are solely generated by RG evolution. Un- 
fortunately we do not have much interesting phenomenol- 
ogy here: the decay K + — > n + vV turns out to have a very 
tiny contribution, orders of magnitude smaller than the 
SM prediction. 
Case 2: AJ 23 ± 

In this case the bounds are controlled by v phe- 
nomenology only: none of the three A^ 23 can be more 
than 2.6 x 10~ 3 at Mq. This yields unobscrvably tiny 
contributions in all rare weak processes. 
Case 3: A' 123 ^ 0, A' 213 , A 313 ^ 

Following the identical iterations as in Case 1, our 
optimum choice for GUT scale inputs come out to be 
A' 123 = 0.00 26, A' 213 = 0.0194, and A 313 = 0.0386. With 
this, all low-energy constraints on product couplings are 
satisfied (including that on K° — K° mixing and Atuk), 
but this generates a large contribution to K + — > 7r + vv\ in 
fact, with such a choice, the bound on |A 313 A 323 |, which 
is 3.5 x 10~ 4 for 300 GeV squarks [15], is violated by a 
factor of 3. Of course, one can scale down the input val- 
ues at Mq, but this shows that RPV couplings that are 
severely restricted by v phenomenology can still generate 
in a correlated way a large amplitude for the K + — > 7r+z/77 
decay, and this channel is worth watching out. However, 
the BRs of various lepton flavor violating r and p decay 
modes turn out to be much smaller than Case 1, of the 
order of 10~ 14 or less, and hence completely uninterest- 
ing. 



V. LEPTON NUMBER VIOLATING DECAYS 
OF SPARTICLES 

A. Benchmark points 

As in the last section we shall consider three choices of 
relatively large trilinear couplings as inputs at the GUT 
scale: 

Case (1) : A^ 13 ^ at M G , 
Case (2) : A^ 23 ^ at M G , 

Case (3) : any combination of three non-zero couplings 
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taken from Case 1 and Case 2. 

In all three cases the chosen trilinear couplings must carry 
different lepton indices i. All the three choices can allow 
direct lepton number violating decays of the sparticles 
(sleptons, squarks and LSP) with appreciable BRs. More 
significantly the collider signals can reveal, both qualita- 
tively and quantitatively, the GUT scale RPV physics 
related to the origin of v masses. 

From the upper bounds on the input couplings the pre- 
cise magnitudes of the input couplings cannot be com- 
puted. However, it has been noted in [31] that the hi- 
erarchy A^ 13 < AJ 2 3,A^33 provides an interesting solution 
of the oscillation data. If we set the magnitudes of the 
three input couplings in Case 1 in the ratio of their upper 
bounds, the desired hierarchy is maintained. For the pur- 
pose of illustration we shall follow this procedure, namely, 
of setting the input RPV couplings at Mq at such values 
so as to reproduce, at M\y, the upper bounds of the re- 
spective couplings computed with a common mass of 100 
GeV for all the sfermions. (Note that they may not be 
the actual upper bounds at the benchmarks chosen by us, 
but one needs a standard set, valid across the different 
benchmark points, to make quantitative comparisons.) 
In case (2) the signal is qualitatively different from case 
(1). Even here, for quantitative estimates we shall set all 
input couplings equal to their upper bounds (i.e., at 100 
GeV). Using A' scaled up by the ratio m^/100 (see Sec- 
tion IV) one obtains larger RPV signals. Our estimates 
are, therefore, conservative. 

Sparticles have RPC decays as well. So the BR of the 
RPV decay of any sparticle depends on the parameters 
of both RPV and RPC sectors. Computing the BRs us- 
ing the above prescription on the input couplings, we 
find that sizable BRs of direct RPV decays of sfermions 
are still allowed. However, for quantitative estimates, we 
show our numbers for some favorable benchmark points, 
for example, the LHC benchmark is at 

too = 200 GeV, toi = 250 GeV, 

A = 0, tan (3= 10, sgn(^i) = -1. (13) 

From radiative electroweak symmetry breaking one gets 
\fi\ = 345.44 GeV. The SUSY soft parameters are at 
Mi = 107.1 GeV and M 2 = 210.4 GeV, and the mass 
spectrum (in GeV) is 

m e - L = m^- L = 270.4, = 220.9, to t ~ 2 = 271.9, m u ~ L = 
576.7, to r = 582.1, m.~ = 526.3, m h ~ = 554.7, m - = 

' »L ' Ol '0 2 ' Xl 

106.4, m -o = 199.7, m- = 350.4, m - = 361.8, m -. = 

' xl ' xl ' xl xt 

200.0, m -. = 365.1, to 5 = 634.7. 

xi y 

We have used three representative benchmark points 
for the Tevatron. 

Choice (1): Signals for gaugino pair production 

too = 200 GeV, mi = 145 GeV, 

A = -540, tan/3 = 11, sgn{p) = -1. (14) 

Radiative electroweak symmetry breaking yields 
\fi\ = —318.9 GeV, and SUSY soft parameters are 



Mi = 61.4 GeV , M 2 = 121.9 GeV. The SUSY mass 
spectrum (in GeV) is 

m e ~ L = m^ L = 229.7, m r ~ = 198.9, mf 2 — 230.7, m,u L — 

381.3, m d ~ L = 389.4, = 313.8, m 6 - = 367.8, m t - = 

153.1, to/- = 406.4, m -„ = 61.3, to - = 117.6, to -+ = 

' * 2 ' x" ' x% xt 

117.6, to ~ + = 337.3, m~ g = 390.4. 

X-2 

This choice is motivated by the fact that the elec- 
troweak gaugino masses are slightly above the kinematic 
reach of LEP and should be observable at the Tevatron. 
Choice (2): Signals from top squark pair produc- 
tion 

Here we use two benchmark points, namely, 

(2a) : too = 140 GeV, toi = 180 GeV, 

A) = -631, tan/3 = 11, sgn(n) = -1, (15) 

and 

(2b) : m = 140 GeV, toi = 180 GeV, 

A Q = -630, tan/3 = 6, sgn(n) = 1. (16) 

For choice (2a) we have = -381.1 GeV, Mi = 76.5 
GeV, and M 2 = 151.32 GeV. The mass spectrum (in 
GeV) is 

m e ~ L = m^ L = 195.1, = 136.1, to t ~ 2 = 200.4, m^ L — 

423.2, m d - L = 430.6, m 6 - = 346.7, m fe - = 407.7, m t - = 

135.5, mr L = 456.9, to = 74.3, to - ' '= 140.9, to + = 

' * 2 ' x? ' xl ' xt 

140.8, m ~+ = 398.0, m~ g = 468.5. 

X2 

The above choice is motivated by the requirement that 
the top squark be the NLSP. For choice (2b), tan/3 is 
smaller so that the RPC loop decay is somewhat sup- 
pressed and the RPV decay with smaller couplings, as in 
Case (2), can compete with the RPC decay [9]. For this 
case = 385.4 GeV, Mi = 76.5 GeV, and M 2 = 151.3 
GeV. The mass spectrum (in GeV) is 
to c ~ l = m^ L = 194.8, TOf x = 152.0, rrif 2 = 198.3, m u ~ L — 
423.2, m d - L = 430.4, to 6 - = 354.1, m 6 - = 410.7, m t - = 

128.4, to.; = 462.2, to = 73.4, to - Q 2 = 138.8, to + = 

2 Xi X2 Xi 

138.4, to - + = 402.9, TOg = 468.7. 

Xi 

All production cross-sections and BRs in this section 
have been calculated using the software CalcHep [33]. 

B. Lightest neutralino decay 

We begin our discussions with RPV decays of the LSP, 
assumed to be Xi ■ These decays govern all RPV signals 
at the colliders except for the ones arising from direct 
RPV decays of sfermions. If there are only A- 33 -type 
couplings, i.e., only those required by v physics, then 
the main LSP decay modes are xi° — * "ibb (assuming 
mt > to^o). As has been discussed in [9] the lepton 
number violating nature of this decay is not obvious due 
to the missing z^s. Moreover, this decay can be faked by 
RPC decays like \2° Xi bb. This is especially so if the 
masses of xi an d X2 do not follow the hypothesis of 
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gaugino mass unification at Mq and the LSP is allowed 
to be much lighter than the x~2°- R should also be noted 
that the above decay may be the main decay channel of 
X2 if the lighter b-squark happens to be significantly 
lighter than the other sfermions due to strong mixing 
effects at large tan (5. 

In our model the LSP decay pattern is totally differ- 
ent because of the presence of the relatively large A- 13 
and/or A^ 23 couplings. The dominant decay modes are 
Xi — * {liu/vid)b (Case 1) or xi —* (hc/vis)b (Case 2). 
In Case 3 some lepton flavors will be accompanied by 
u/d jets while the others will come in association with 
c/s jets. Thus the nature of the heavy flavor jets accom- 
panying the charged lepton in the final state would be 
the qualitative feature that can distinguish one choice of 
GUT scale physics from another. 
Case 1: ? 

Here our choices are A' 113 = 0.0064, A' 213 = 0.0194, 
and A 313 = 0.0386 at Mq (see Section IV), which we use 
for the purpose of illustration in our numerical compu- 
tations. At the weak scale the magnitudes increase by 
about a factor of 3. For our model of quark mixing (Sce- 
nario I) the physical couplings for the decays involving vs 
are the same as the ones in the flavor basis. For the de- 
cays into charged leptons the couplings in the two bases 
are practically identical as they are related by diagonal 
elements of the CKM matrix. 

These couplings generate several novel channels for xi° 
decay, tabulated in Table II. We have also presented in 
this table the partial widths and the BRs of different 
modes. The BRs for other values of the couplings can be 
computed by simple scaling of the partial widths. It is to 
be noted that the lepton number violating nature of the 
underlying interaction is obvious from the decays involv- 
ing charged leptons. Since, Xi 1S a Majorana particle it 
can decay into leptons of both charges with equal BRs. 
Thus final states with like sign dileptons or dileptons of 
different flavors and no missing energy would be clean 
signals of sparticle pair production in this model. More- 
over, the relative abundance of final states with different 
leptonic compositions will be a strong indication of the 
magnitudes of underlying input couplings at Mq- 

Finally from the total width computable from Table 
II, we find that the LSP is going to decay inside the 
detector. Note that it is only a matter of proper scaling 
to get the partial widths with a different set of input 
couplings. However, for larger input couplings, the Xi 
lifetiem will be correspondingly shorter and the displaced 
vertex may be missed. 

It is interesting to note that in addition to the highly 
suppressed A^ 33 couplings, required for ^-physics, this 
choice of input couplings also induces at M\y A'i23 cou- 
plings with reasonable magnitudes via CKM rotation, as 
shown in Table II. The decays triggered by these induced 
couplings into final states involving charged leptons and 
charm have small BRs, but they may still be observable 
in the clean environment of the ILC. 
Case 2: AJ 23 ± 



X 1 v-l CLy o 




l_y CL V 


Pq r+ i a 1 
_L dL Lid! 


BR 
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at M w 
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\i 

A 113 


n no 


e ub 


1 w 1 r\— 9 
l.Z X 1U 


0.0 X 1U 






i/ e db 


1.4 X 10~ 9 


1.0 x 10~ 2 


A213 


0.064 


pTub 


1.3 x 10 -8 


9.0 x 10 -2 






TT^db 


1.4 x 10 -8 


1.0 x 10 _1 


A313 


0.12 


r~ub 


5.4 x 10 -8 


3.8 x 10" 1 






u^db 


5.3 x 10" 8 


3.8 x 10 _1 


A123 


0.0045 


e~ cb 


6.3 x 10- 11 


4.5 x 10~ 4 


V23 


0.014 


[i~ cb 


6.1 x 10~ 10 


4.3 x 10~ 3 


-^323 


0.03 


t~ cb 


3.4 x 10~ 9 


2.4 x 10~ 2 



TABLE II: BRs of LSP decays at LHC, for Case (1) parameter 
set. The CP conjugate channels are implied. The last three 
couplings are generated by CKM rotation and can produce 
only charged leptons in the final state. 



Xi u decays 


Value 


Decay 


Partial 


BR 


through 


at Mw 


channel 


width (GeV) 


at LHC 


^'123 


0.0082 


e~ cb 


2.1 x 10~ 10 


1.5 x 10 _1 






TT^sb 


2.3 x 10" 10 


1.6 x 10" 1 


A223 


0.0082 


pT cb 


2.1 x 10 -10 


1.5 x 10 _1 






TT^sb 


2.3 x 10" 10 


1.6 x 10" 1 


A323 


0.0082 


t~ cb 


2.5 x 10" 10 


1.8 x 10" 1 






Ts^sb 


2.5 x 10" 10 


1.8 x 10" 1 


A113 


0.0018 


e~ub 


1.0 x 10~ n 


7.1 x 10~ 3 


Vl3 


0.0018 


pT ub 


1.0 x 10~ n 


7.1 x 10~ 3 


A313 


0.0017 


T~ub 


1.1 x 10 -11 


7.7 x 10 -3 



TABLE III: Same as Table II, but for Case (2) parameter set. 



Xi u decays 
through 


Value 

at M w 


Decay 
channel 


Partial 
width (GeV) 


BR at 
Tevatron 


A113 


0.02 
(0.02) 


e~ub 
T^db 


2.2 x 10~ 10 
(6.3 x 10- 10 ) 
2.2 x 10~ 10 
(1.1 x 10" 9 ) 


1.0 x 10~ 2 
(5.6 x 10~ 2 ) 

1.3 x 10~ 2 
(9.8 x 10~ 2 ) 


-^213 


0.064 
(0.027) 


\i~ub 
T^ldb 


1.9 x 10" 9 
(1.2 x 10~ 9 ) 

2.2 x 10~ 9 
(2.0 x 10~ 9 ) 


9.1 x 10" 2 
(1.0 x 10 _1 ) 

1.1 x 10 _1 
(1.7 x 10" 1 ) 


A313 


0.12 
(0.03) 


T~ub 
TT^db 


8.0 x 10" 9 
(3.2 x 10~ 9 ) 

7.9 x 10~ 9 
(2.9 x 10" 9 ) 


3.9 x 10 _1 
(2.9 x 10" 1 ) 

3.8 x 10" 1 
(2.6 x 10 _1 ) 


V23 


0.0045 
(0.0046) 


e~ cb 


1.1 x 10~ n 
(3.3 x 10 -11 ) 


5.3 x 10~ 4 
(3.0 x 10 -3 ) 


V23 


0.014 
(0.0061) 


u~ cb 


1.1 x 10" 10 
(5.9 x 10 -11 ) 


5.0 x 10 -3 
(5.2 x 10 -3 ) 


Vi23 


0.03 
(0.0065) 


r~cb 


5.2 x 10 -10 
(1.7 x 10 -10 ) 


2.5 x 10~ 2 
(1.5 x 10~ 2 ) 



TABLE IV: Same as Table II, for Case (1) and Tevatron 
benchmark point 1 (point 2(a)). 
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Xi u decays 
through 


Value 
at Mw 


Decay 
channel 


Partial 
width (GeV) 


BR at 

_1_ Mi V CIU1 \JIL 


■^123 


0.0082 


e~cb 
TT^sb 


3.6 x lO^ 11 
(9.1 x 10~ n ) 
5.0 x 10~ n 

(Z.2 X 10 J 


1.4 x 10" 1 
(8.8 x 10~ 2 ) 

2.0 x 10 _1 
(z.l X 10 J 


^223 


0.0082 


[i~ cb 
T^sb 


3.6 x 10" 11 
(9.1 x 10" 11 ) 
5.0 x 10~ n 

1 Z.Z A 1U 1 


1.4 x 10 _1 
(8.8 x 10~ 2 ) 

2.0 x 10 _1 
(o i v i n — ^\ 

1 Z.l A 1U 1 


^323 


0.0082 


r~cb 
T^sb 


3.5 x 10 _ii 
(1.3 x 10~ 10 ) 

4.2 x 10~ n 
(2.6 x 10" 10 ) 


1.4 x 10" 1 
(1.2 x 10 _1 ) 

1.7 x 10 _1 
(2.5 x 10 _1 ) 


-^113 


0.0018 


e~ub 


1.7 x 10~ 12 
(4.4 x 10~ 12 ) 


6.8 x 10 -3 
(a 3 x 10 -3 ) 


^213 


0.0018 


fj,~ub 


1.7 x 10 -12 
(4.4 x 10" 12 ) 


6.8 x 10 -3 
(4.3 x 10 -3 ) 


■^313 


0.0018 


T~ub 


1.6 x 10" 12 
(6.2 x 10~ 12 ) 


6.4 x 10" 3 
(6.0 x 10~ 3 ) 



TABLE V: Same as Table II, for Case (2) and Tevatron bench- 
mark point 1 (point (2b)). 



In this case none of the three input couplings X' i23 can 
be larger than 2.6 x 1CT 3 at M G (see Section IV). The 
allowed LSP decay modes, shown in Table (III), are char- 
acterized by the presence of charm and strange particles 
in the final state. These jets will be isolated from the b jet 
and hopefully can be tagged. Moreover the strangeness 
quantum number of the strange particle directly emit- 
ted by the Xi will be opposite to that emitted by the b 
(unless that hadronizes into a neutral B meson and os- 
cillates). It is encouraging to note that the possibility of 
c-jet tagging and reconstruction of c-flavored hadrons are 
being discussed vigorously [34]. 

This particular set of input couplings will induce A- 13 
with much reduced strengths compared to the input val- 
ues in case 1. Thus the dominant (rare) decay modes 
in this case would be the rare (dominant) decay modes 
of case 1. The magnitudes of the induced couplings are 
given in Table III. 
Case 3 

Here the three input couplings could be various combi- 
nations of the six couplings presented in Case 1 and Case 
2 where each input coupling should bear a different lep- 
ton index. The decay modes will now be combinations of 
the ones presented in Tables II and III. Using the partial 
widths presented in these tables one can easily compute 
the desired BRs. Again the decay modes of the LSP and 
the relative population of different final states may lead 
to the underlying model. 

In Tables IV and V, we have shown the LSP decay 
modes for the Tevatron benchmark points. Note that for 
Case 2, we have chosen benchmark 2(b), the reason of 
which has been explained earlier. 



C. Direct RPV decays of the sleptons 

With our choices of A' 113 and A 213 in Case 1 the direct 
RPV decays of the left selectron ez, and the left smuon 
fiij have strongly suppressed BRs (see Table VI). Such 
decays may therefore occur only as rare modes. In sharp 
contrast the heavier r slepton mass eigenstate, which is 
dominantly a left slepton, may have sizable BRs for di- 
rect lepton number violating channels. However, if we 
scale the input A' couplings by m^/100 for the LHC 
benchmark point, even and fiL will have RPV BRs 
competitive with the RPC channels. In case 2 the direct 
RPV decays of all sleptons will only occur as rare modes. 



D. Direct RPV decays of the squarks 

There is a small but non-negligible parameter space 
within the mSUGRA framework where gluinos are heav- 
ier than all squarks. Our choice of mSUGRA parameters 
belongs to this parameter space. However, this condi- 
tion may be satisfied by squarks belonging to the third 
generation in a larger region of the parameter space. In 
the mSUGRA model the b squark mass eigenstates can 
be lighter than the gluinos at large or even intermediate 
values of tan/3 due to mixing effects. The more interest- 
ing case of the top squark will be discussed seperately. 

Squarks lighter than the gluinos are more common if 
non-universal squark and/or gluino masses [30, 35] moti- 
vated by various intricacies of physics at Mq are allowed. 
Detailed phenomenology of RPC models with squarks 
lighter than the gluinos due to SUSY breaking SO(10) 
D-terms has been discussed in [36] . The RPV signals can 
be obtained from these by simply adding the LSP decay. 
Since the RG equations of the RPV A'-type couplings 
are almost independent of the RPC parameters, our es- 
timates of the magnitudes of these couplings remain by 
and large unaltered even in models beyond mSUGRA. 

In the LHC benchmark point and in Case (1), squarks 
dominantly decay into electroweak gauginos. We present 
in Table VII the BRs of RPC decays, as well as those 
of RPV decays through A- 13 couplings of ul and cIl- 
Only the decays into r or v T may have a few percent 
BR while the other modes are suppressed. Table VIII 
shows the channels for b squark decay. For the lighter 
mass eigenstate b\ RPV modes have highly suppressed 
BRs, while the heavier mass eigenstate 62 can decay into 
several RPV channels each having a few percent BR. It 
should be noted that the BRs of the RPV decays of b 2 
will increase dramatically if < m t + . 

E. Direct RPV decays of the lighter top squark 

The lighter top squark mass eigenstate t\ can natu- 
rally have a mass much smaller than that of the other 
squarks. This may happen due to two reasons. First, 
even if all squarks have the same mass mo at Mq, the 
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Tvne of 


RPC Decav 


Partial 


Branching 


RPV Decay 


Partial 


Branching 


Slepton 


Channel 


Width (GeV) 


Ratio 


Channel 


Width (GeV) 


Ratio 


e~L 


Xi e 


z.o X 1U 


Z.O X 1U 










A.2 c 


2.5 x 10 _1 


2.7 x 10 _1 


ub 


6.6 x 1CT 3 


7.2 x 1CT 3 




Yl ~~ 


4.3 x 10 _1 


4.7 x 10" 1 


cb 


3.3 x 1CT 4 


3.6 x KT 4 




Xi A 4 


o o w in — 1 
1.6 X 1U 


o /i w in — 1 
2.4 X 10 










A^ A* 


2.5 x 1CT 1 


2.5 x 10 _1 


ub 


6.7 x 1CT 2 


6.8 x 1CT 2 




A 1 M 


4.3 x 10 _1 


4.4 x 1CT 1 


cb 


3.2 x 10~ 3 


3.3 x 10~ 3 


fi 


Xl T 


r i i n~ i 
0.1 X 1U 


9. / X 1U 


ub 


i i i n~2 
1.1 X 1U 


i 1 in — 2 

i. ( x iu 




X2 ' 


2.6 x 1CP 3 


4.1 x 10~ 3 


cb 


6.6 x 10~ 4 


1.1 x 10 -3 




Xi"*^ 


4.6 x 10~ 3 


7.3 x 10" 3 










xA 


2.7 x lfT 1 


2.4 x 10 _i 


ub 


2.2 x KT 1 


2.0 x 10 _I 




X2°r 


2.2 x 10 _1 


2.0 x 10 _1 


cb 


1.4 x 1(T 2 


1.3 x 10~ 2 




Xi "t 


3.7 x KT 1 


3.4 x 10" 1 









TABLE VI: BRs of RPC and RPV decay channels of sleptons in Case (1) with the LHC benchmark point. 



Type of 


RPC Decay 


Partial 


Branching 


RPV Decay 


Partial 


Branching 


Squark 


Channel 


Width (GeV) 


Ratio 


Channel 


Width (GeV) 


Ratio 




Xi U « 


6.3 x 10~ 2 


1.1 x 10~ 2 


e + b 


4.6 x 10~ 3 


8.1 x 10~ 4 




X2°U 


1.7 


3.1 x 10 _1 








U~L 


X3% 


9.0 x 10~ 3 


1.6 x 10~ 3 


fi + b 


4.7 x 10~ 2 


8.3 x 10~ 3 




X~4°« 


8.1 x 10~ 2 


1.4 x 10~ 2 










Xi + d 


3.5 


6.2 x 10 _1 


r+b 


1.6 x 10 _1 


2.9 x 10~ 2 




X2 + rf 


7.3 x 10~ 2 


1.3 x 10" 2 










Xi"d 


8.4 x 10~ 2 


1.5 x 10~ 2 


VZb 


4.0 x 10~ 3 


7.2 x 10~ 4 




X2°d 


1.7 


3.1 x 10 _1 










X3°d 


1.5 x 10~ 2 


2.7 x 10~ 3 




4.1 x 10~ 2 


7.5 x 10~ 3 




X4°d 


1.1 x 10 _1 


2.0 x 10~ 2 










Xi w 


3.1 


5.6 x 10" 1 


IT^b 


1.4 x 10 _1 


2.6 x 10~ 2 






2.9 x 10 _1 


5.5 x 10~ 2 









TABLE VII: BRs of RPC and RPV decay channels of up and down squarks in Case (1) with the LHC benchmark point. 



mass parameters of it and tn at the weak scale may 
be significantly smaller due to the influence of the large 
top quark Yukawa coupling in the RG equations. The 
mass of the lighter eigenstate may be further suppressed 
due to mixing effects in the t^-tn mass matrix, with 
large off-diagonal entries, i It is quite conceivable that t\ 
happens to be the next lightest supersymmetric particle 
(NLSP). If in addition the conditions < m t + m^o 
and mi t < mi, + miy + m^o are satisfied, then the only 
possible RPC decay modes are 
(i) the loop decay [37]: t\ — > cxi° 

or (ii) the four body decay [38]: t\ — > bff xi°, where 
/ and /' are light fermions. Both the decays occur in 
higher order of perturbation theory and have naturally 
suppressed widths. If the couplings of the type A- 3 j are 
nonzero and of the order of 1CP 1 or 10~ 2 , then the RPV 
two body decay 

ti - l+d] (17) 

overwhelm the RPC decays and occur with 100% BRs. 
However, If the underlying trilinear coupling is order 



10~ 3 or smaller then the BR of each competing mode 
may have a sizable magnitude. In fact it was shown in 
[8] that the data from Tevatron Run I or the preliminary 
data from Run II are already sensitive to A^ 3j ~ 10~ 3 - 
10~ 4 , although for rather small m^. It was shown in [9] 
that a larger range of can be probed at Run II with 
integrated luminosity ~ 2 fb _1 . 

If A^ 33 , i.e., the couplings required by the v sector, are 
the only ones with appreciable magnituds then the main 
decay channels of t\ would be 

*i - in. (is) 

In contrast our model allows other possibilities, particu- 
larly in Scenario II (Section VI), depending on the un- 
derlying model of CKM mixing and the input RPV cou- 
plings. First consider Scenario I. Here the three input 
couplings in the flavor basis at Mq are possible combi- 
nations of A^ 13 and Aj 2 3 ■ These input couplinngs can 
induce A - 33 at the weak scale via the RG evolution as we 
have already discussed. However, the input couplings, 
depending on their magnitudes, may induce couplings of 
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lype or 


i v i l uecay 


Partial 


Branching 


txr v uecay 


Partial 


Branching 


Scjuark 


Chann6l 


VV lQlll ^Vjr6V ) 




Channel 


vv iaxn ^Vje v ) 


ixario 




- (It 

Xi b 


1.7 x 10~ x 


4.8 x 10~ 2 


e~u 


5.6 x 10" 4 


1.6 x 10~ 4 




X2°b 


1.5 


4.1 x 10 _1 


flU 


5.8 x 10~ 3 


1.6 x 10~ 3 


h 


Xs°b 


8.8 x 1(T 2 


2.5 x 10~ 2 


TU 


2.0 x 10~ 2 


5.7 x 10~ 3 




X4°b 


8.2 x 1(T 2 


2.3 x 10~ 2 


VZd 


4.8 x 10~ 4 


1.4 x 10~ 4 




Xi t 


1.7 


4.8 x 10 _1 


T^d 


4.9 x 10~ 3 


1.4 x 10~ 3 




X~2 t 


does not 


occur 


T\d 

ec 

flC 
TC 


1.7 x 10~ 2 

2.8 x 10~ 5 
2.8 x 10~ 4 
1.3 x 10~ 3 


4.8 x 10~ 3 

7.9 x 10~ 6 
7.9 x 10~ 5 
3.6 x 10~ 4 




~ 1 1 1 

Xi b 


2.0 x 1CT 1 


1.3 x 10 _i 


eu 


3.8 x 10~ 3 


2.4 x 10~' 1 




X2°b 


7.4 x 1(T 2 


4.7 x 10~ 2 


flU 


3.9 x 10~ 2 


2.5 x 10 -2 


b-2 


Xs°b 


1.8 x 10 _1 


1.1 x 10 _1 


TU 


1.4 x 10 _1 


8.8 x 10~ 2 




X4°b 


1.9 x KT 1 


1.2 x 10 _1 


TJld 


3.3 x 10~ 3 


2.1 x 10~ 3 




Xi t 


9.0 x 10~ 2 


5.8 x 10~ 2 


TJJId 


3.4 x 10~ 2 


2.2 x 10~ 2 




X.2 t 


4.8 x 10 _1 


3.1 x 10 _1 


Twd 

ec 

)1C 
TC 


1.2 x 10 _1 
1.9 x 10~ 4 
1.9 x 10~ 3 
8.6 x 10~ 3 


7.6 x 10~ 2 
1.2 x 10~ 4 
1.2 x 10~ 3 
5.5 x 10" 3 



TABLE VIII: BRs of RPC and RPV decay channels of bottom squarks in Case (1) with the LHC benchmark point. 



the type A- 33 involving charged lepton interactions via 
the CKM rotation which could be much larger than the 
same couplings relevant for the v sector via RG evolution 
only. This is a more likely situation for input couplings 
with possibly larger magnitudes as in Case 1. 

These A- 33 couplings trigger top squark decay. Maxi- 
mum possible values of these couplings at the weak scale 
are: 



A' 133 = 8.0xl(T 4 , A2 33 = 2.4xl(T 3 , A 333 = 4.6 x 10~ 3 . 

( 19 ) 

The RPV decay modes will still be t\ — > l+b. The com- 
plete dominance of RPV decays, however, would strongly 
indicate the relatively large input couplings at the GUT 
scale. 

For the purpose of illustrating the competition between 
RPC and RPV decay modes, we prefer to use somewhat 
scaled down values of the RPV couplings of Case (1) 
[otherwise the RPV decay modes will be overwhelmingly 
large] . The input values as well as the BRs are presented 
in Table IX. For Case (2), on the other hand, the up- 
per limits of the input couplings A - 23 are much smaller 
in magnitude (~ 10~ 3 ), and the CKM rotated A i33 cou- 
plings generated by them would be naturally ~ 1CP 4 or 
smaller. Thus the possibility of competition among RPC 
and RPV is better in this case, see Table X. We use 
the Tevatron benchmark values, but let us emphasize 
that the stop signals will be much more prominent at 
the LHC. 

In scenario II the situation is even more intriguing. For 
a discussion we refer the reader to Section VI. 



Type of 


Decay 


RPV coupling 


Branching 


Squark 


Channel 


and strength 


Ratio 




ti -» Xi'c 




15.5% 




ii -> e+b 


A' 133 : 8.0 x 10- 4 


18.5% 


tl 


ti ->n + b 


A 233 : 1.03 x 10~ 3 


30.7% 




tl -> T+b 


A 333 : 1.1 x 10" 3 


35.0% 



TABLE IX: The lighter top squark BR into RPC and RPV 
channels. The input is as in Case (1), with couplings at the 
mass basis and at Mw ■ The benchmark point is that of Teva- 
tron (2a). 



Type of 


Decay 


RPV coupling 


Branching 


Squark 


Channel 


and strength 


Ratio 




ti Xi'c 




34.2% 




ti -^e+b 


A'133 : 1.3 x 10~ 4 


21.8% 


ti 


ti -> [i+b 


A~2 33 : 1.3 x 10~ 4 


21.8% 




tl -> T+b 


A 333 : 1.3 x 10" 4 


21.8% 



TABLE X: The same as Table IX, but with Case (2) input 
and Tevatron (2b) benchmark. 



F. Some interesting signals with sizable number of 
events at the Tevatron and the LHC 

Our aim is to propose some signals, in addition to 
LSP decays, which indicate the direct RPV decay of 
at least one sfermion. A signal arising from sfermion- 
antisfcrmion pair production followed by direct RPV de- 
cays of both is not suitable for this purpose because of 
obvious backgrounds. Slepton pair production, for ex- 
ample, would lead to final states with jets only (see Ta- 
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Channel 
number 


Decay modes of 
stop pair 


No. of 
events 


Channel 
number 


Decay modes of 
stop pair 


No. of 
events 


i. 


ti -> r+b 
ti — ► e o 


3403 


2. 


ti r+b 
ti -^(ib 


5647 


3. 


ti ^i+b 
ti — ► e o 


2985 


4. 


ti r+b 
ti — > xi c — * e be u 


80 


5. 


ti -> r+b 

tl — > X~1°C — * T+bc U 


413 


6. 


ti e+b 
ti — > Xi°c — > j-i+ be u 


75 


7. 


tl /i+b 

ti — » xi c — > /i+bc w 


125 


8. 


ti -> ^+b 
ti* — > xi c — ► r+bc u 


363 


9. 


ti -» e+b 
ti* — > Xi°c — » r+bc u 


218 


10. 


ti -f r+b 
ti* — > Xi°c — * /i+bc u 


142 



TABLE XI: Number of events coming from stop-antistop pair production at Tevatron for final states involving dileptons with 
(i) like sign and same flavor; (ii) like sign and different flavor; and (iii) unlike sign and different flavor. The dilepton signal 
is accompanied by two jets. The charge conjugated states are also included. The input parameters are of Case (1), Tevatron 
benchmark (2a), with m tl = 135.5 GeV. The stop pair production cross-section is 2.92 pb. We use the data of table IX and 
IV. 



Channel 
number 


Decay modes of 
stop pair 


No. of 
events 


Channel 
number 


Decay modes of 
stop pair 


No. of 
events 


1. 


ti -> T+b 

ti — > e~b 


3482 


2. 


ti -> T+b 
ti — ► (i~ b 


3482 


3. 


ti -» n+b 

ti — ¥ e~b 


3482 


4. 


ti -» e+b 
ti* — ► x"i°c — ► /i+bc c 


240 


5. 


tl -> T+b 

ti* — > Xi"c — ► T+bc c 


328 


6. 


tl — > Xi° c — * i~+bc c 
ti* — > x"i°c — ► r+bc c 


31 


7. 


ti -> /k+6 
ti* — > Xi°c — ► /i+bc c 


240 


8. 


ti — > xV'c ~* /i+bc c 
ti* — ► x~i°c — ► /i+bc c 


17 


9. 


ti -> e+b 
ti* — > xi°c — > e+bc c 


240 


10. 


ti — > Xi U c — * e+bc c 
ti — > xi c — > e+bc c 


17 


11. 


ti -> /k+6 
ti* — > x"i°c — ► r+bc c 


328 


12. 


ti -> e+b 
ti — > xi c — ► r+bc c 


328 


13. 


ti -f e+6 
ti* — > x"i°c — ► r+bc c 


240 


14. 


ti -> r+b 
ti* — ► Xi°c — * /t+bc c 


240 



TABLE XII: Same as Table XI, but for Case (2) and Tevatron benchmark (2b), where rh tl = 128.4 GeV, and the stop pair 
production cross-section is 4.07 pb. We use the data of table X and V. 



ble. VI) which will suffer from a huge QCD background. 
Squark-antisquark pairs decaying into isolated opposite 
sign dileptons (OSDs) of the same flavor and accompa- 
nied by 66 jets may suffer from a large Drell-Yan back- 
ground accompanied by QCD jets . However, final states 
with isolated like sign dilepton (LSD) pairs or dileptons 
carrying different flavors are expected to provide cleaner 
signals. 

We discuss below some interesting signals arising from 
sparticle pair production and their subsequent decay. To 
estimate the number of events, we use an integrated lu- 
minosity of 10 5 pb- 1 at LHC and 9000 pb" 1 at Tevatron. 

1. Signals from sfermion pair production 

First we consider the dilepton plus dijet signals from 
direct RPV decays of the lighter top squark pair pro- 



duced at the Tevatron. Even for the apparently unfa- 
vorable signal involving opposite sign dileptons (OSD) of 
the same flavor, the backgrounds can be suppressed by 
suitable cuts [8, 9]. We do not analyze this signal in this 
work. Rather, we analyze only those channels which re- 
sult in either OSD of different flavor or LSD in the final 
state. We present in Tables XI and XII the number of 
events using the BRs in Tables II, III, IX, X. The num- 
ber of events in several cases are quite encouraging. The 
reconstruction of the masses of the two top squarks when 
both decays directly to RPV channels was discussed in 
[9]. In some of the cases discussed, we need the recon- 
struction of the invariant masses of two sets: one with a 
lepton and a 6-jet, and one with all other visible particles 
in the final state, such that the two reconstructed invari- 
ant masses agree with each other within certain tolerance 
limit. This establishes our signal. 

Next we focus our attention on the final states 
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ij ± Zj ± &6(66)+X, where X is the unobserved junk at LHC 
arising from initial states not involving top squarks. One 
of the isolated leptons in the like sign dilepton (LSD) 
pair comes from the RPC decay of a sfermion followed 
by the RPV decay of the produced LSP, which is a Majo- 
rana fcrmion and decays into leptons of both flavors with 
equal BRs. The other one comes from the direct RPV 
decay of the second member of the produced fermion- 
antisfcrmion pair. Such final states may arise from sev- 
eral initial states listed in Tables XIII and XIV. The 
background to this signal is expected to be small due to 
the presence of the isolated LSDs without any accompa- 
nying missing energy. In principle the presence of the 
bb or bb pair in the final state may also help to reduce 
the QCD background which always involve bb. This is 
possible if lepton tagging and/or kaon tagging can dis- 
tinguish between b and b with reasonable efficiency. One 
can hope that this will be achieved since the study of CP 
violation in B-decays is an important component of LHC 
physics. The direct lepton number violating decay of one 
of the sfermions can be identified by reconstructing the 
invariant mass as in the case of the decaying top squarks. 

From Table. VII it is clear that the squarks may domi- 
nantly decay into rs via the RPV channel while decays to 
e and fi may be heavily suppressed. However, the BR of 
the muonic channel is still of the order of a few percent. 
Combining these informations one can estimate that a 
significant number of events of r — fi type (of all charge 
combinations) can be seen. We have also presented in 
Table XIII (see last two rows) r ± /i ± 66 events which are 
also expected to have small backgrounds. The signal size 
is already encouraging. 

The production cross-section of slepton pairs is small 
at LHC. In spite of this, the signal may be important 
if the signal from the squarks are depleted due to the 
presence of gluinos lighter than the squarks. The LSD 
signal from r slepton pair production has apparently a 
small size, but this may be enhanced using larger values 
of A' which are still allowed (so this is in some sense a 
conservative estimate). 

2. Signals from gaugino pair production 

The T + r + bb signals from gaugino pair production at 
LHC and Tevatron are presented in Table XV for Case(l) 
and in Table XVI for Case (2) input values. Sizable 
number of events are expected from X2°Xi + an d Xi + Xi~ 
pairs involving LSDs of a particular flavor. 

We conclude this section by reiterating that all the 
input A's at Mg can in principle be determined by a fit to 
the oscillation data if at least a partial information about 
the RPC sector is available from collider experiments. 
This in turn will firmly predict the particle content of 
different final states as well as the relative abundance of 
various final states. The observation of these exciting 
signals will reveal the GUT scale physics lying at the 
origin of v masses and mixing angles. 



VI. MIXING IN THE DOWN QUARK SECTOR: 
PHENOMENOLOGY 

Let us now discuss how far the neutrino data can be 
useful in constraining the RPV models in Scenario II 
with mixing in the down-quark sector, V = Dt, = Dr, 
Ul = Ur = 1. In this case the off-diagonal elements 
of the Yukawa coupling matrix Yd are much larger than 
those in Scenario I (mixing restricted to the up quark sec- 
tor only) . As a result any input A'-type couplings would 
lead to rather large values of Ki and A- 33 via RG evolu- 
tion, in conflict with the upper bounds in [18] unless the 
magnitudes of the input couplings are too small to be 
of any phenomenological interest. This conclusion agrees 
with that of [20] although the latter was arrived at in a 
model with a single neutrino mass. 

However, this argument hinges crucially on the RG 
evolution from the GUT scale to the weak scale in the 
mSUGRA model. In practice, the origin of SUSY break- 
ing is not known. It is quite possible that SUSY is broken 
at a much lower scale (as, e.g., in the Gauge-mediated 
SUSY breaking models). For this type of models, not 
only the boundary conditions of the RG will be different 
but as a whole the evolution should have a less important 
role to play. 

In view of the above uncertainties we take a more phe- 
nomenological point of view and assume that only a few 
relatively large A' type couplings in the weak basis, not 
directly related to v physics, and bilinear parameters Ki 
are generated at the weak scale by some new physics at 
the high scale, the nature of which is unknown. We will 
see whether it is possible to generate the suppressed A- 33 
type couplings relevant for the v sector via CKM rota- 
tion while the input couplings at the weak scale can still 
be sufficiently large to trigger interesting non-neutrino 
phenomenology 

In this case, the rotation to the mass eigenbasis is dif- 
ferent from that in Scenario I. After rotation eq. (1) reads 

Wf = \ mv N l B m B c p + % jp EiUjD;, (20) 
where the fields are in the mass basis, and 

v v v v* 

^imp ~ /y ijk v ]m^kp^ 
\jp = KjkYkp- (21) 

Thus in contrast to Scenario I new couplings involving 
both neutral and charged lepton supcrficlds can be in- 
duced by the input couplings via CKM rotations. The 
next task is to find out the possible minimal sets of input 
couplings at the weak scale which can serve our purpose. 

Clearly, one can take X' i33 ^ as the input, each one 
with upper limits of 1.5 x 10~ 4 approximately [18]. How- 
ever, the induced couplings in this case will be so small 
that no interesting phenomenology outside the neutrino 
sector apart from LSP decay [32] and top squark decay 
[9] is viable. 

Before we investigate other options, let us note that 
K° -K° mixing data predicts |A^ 12 A^ 21 | < 10~ 9 [13, 15] 
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Channel 
number 


Squark 
pair 


a (pb) 
at LHC 


Decay mode of 
the squarks 


No of 
events 


1. 


— = — — 1 

&2&2 


0.54 


b-2 — * T~U 

62 — * &Xi° — ► T~ub & 


86 


2. 


fe 2 62* 




&2 — > 

bV — > 6yo° — > 7fvi — > r~ub blf 


85 


3. 


bi &r 


0.77 


bi — ► 

hi * — > 6vo° — > T f yi — > r~it6 6 7 f 

u l ' U A2 J J A.1 ' UjU u J J 


69 


4. 


ULUL* 


0.61 


UL* — * T~b 

ul > «X2 — *• / Jxi ^ r u0 J J 


208 


5. 


— = — — 

6262 




&2 — » M + W 
62 — > T It 


237 


6. 


ULUL* 




— > r~b 
u~l — > /i + 6 


32 


7. 


UlUl 


1.02 


iti, — > r + 6 
wi H + b 


49 



TABLE XIII: Number of events arising from squark-antisquark pair production at LHC, with Case (1) input. Charge conjugate 
channels are included. 



Channel 
number 


Slepton 
pair 


a (pb) 
at LHC 


Decay mode of 
the sleptons 


No of 
events 


1. 


T2V2 


5.2 x 10~ 3 


T2* —* Ub 

T2 — ► T~xi° — > t~ t~m& 6 


10 


2. 


T2*T 2 




T2* — > Ub 

t 2 -*■ r"X2° -» 7 /Xi° -» r^itb 6 // 


8 



TABLE XIV: Same as Table XIII, for stau pair production. 



(we drop the tilde from now, but these couplings are in 
the mass basis). Thus, one should not start with nonzero 
Km K22' K12 or K21 ( m the weak basis) as inputs. From 
cq. (21), it is clear that one must start with such small 
values as to be devoid of all phenomenological interests. 
The limiting value of all these four sets is 1.4 x 1CP 4 at 
the weak scale [18]. This translates to |A- 33 | < 2.8 x 
1CT 9 , 2.5 x 1CT 7 , and 2.6 x 1CT 8 (for both X' ll2 and \' l21 
sets) respectively. Considering the symmetric nature of 
the magnitude of the CKM elements, one can start with 
either (i) A^ 13 , i31 j or (ii) A^ 23 , i32 s to be nonzero. At first, 
the couplings are set to be nonzero at the weak scale with 
such a value as to produce the limiting value of 1.5 x 1CP 4 
for A^ 33 , and then rescaled, if necessary, so as to conform 
with the Attik constraint (this is a similar iteration as 
was done in Section II). 

Case (i): The maximum value of A- 23 ^ 32 ^ to start with 
is 3.6 x 10~ 3 . After proper rotation, one obtains the 
couplings at the mass basis and finds |A- 12 A- 21 | < 4.5 x 
1CP 10 , so this choice is safe from the K° — K° constraint. 

Case (ii): The neutrino data constrains A^ 13 , i31 N to be 
less than 0.052, but one needs a scaling down by a factor 
of 9.7 to be in conformity with the K° — K° constraints. 
For all these cases every other constraint on product cou- 
plings is satisfied. 

However, with such small couplings, one does not hope 



to see an indirect non-neutrino signal of RPV. We men- 
tion again that for indirect signals of RPV, the individual 
couplings should be of the order of 10~ 2 and the product 
couplings of the order of 10 -3 -10~ 4 . Such possibilities 
do not exist for the option we have just discussed; but 
again, this is a model-dependent statement and should 
not hold for other models of CKM mixing. 

Thus we conclude that irrespective of the details of 
the high scale physics it is impossible to generate the 
A^ 33 couplings of the benchmark scenario from relatively 
large, phcnomcnologically interesting input couplings at 
the weak scale and in the weak basis if (i) quark mixing 
is primarily restricted to the down quark sector, and (ii) 
the Yukawa matrix is symmetric. 

Still, these couplings can induce novel channels of de- 
cay for the lighter stop quark and the lightest neutralino 
which are quite distinct for those in Scenario I. Two possi- 
ble phenomenological inputs couplings in the flavor basis 
at the weak scale are A^ 31 or A^ 32 . Since these couplings 
themselves can trigger top squark decays, the physical 
couplings in the mass basis does not suffer any further 
suppression due to CKM rotation. The absence of bot- 
tom jets in the final state distinguishes the top squark 
decay from that in Scenario I. These signals have also 
been discussed in refs [8, 9]. 
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Channel 
number 


Gaugino 
pair 


cr (pb) at LHC 
(Tevatron) 


Decay channel of 
the gauginos 


No of events 
at LHC (Tevatron) 


1. 


Xi U Xi° 


9.3 x 1(T 3 


Xi° — » t~u b 
Xi° —> t u b 


62 


2. 


X2 U X 2 ° 


2.1 x 1(T 2 


X2° - / /Xi u - t-u6 / / 

X2 — » / /Xl ^ r uo J J 


130 


3. 


X2°X1 + 


5.3 x 10" 1 

1 _■ . O /\ -L U J 


XVW fxi" ^ r-vb_fj 

V -i V f /n "V 1 ^ > "T H hi ~f f'l 


3350 


4. 


Xi + Xi" 


4.7 x KT 1 
( 3.2 x 10" 1 ) 


Xi + ^/i Mi^r-ubf, h 
Xi /i / 2 Xi° * ^~ "6 A / 2 


3100 
(214) 


5. 


X~2 + X~2~ 


2.3 x 1(T 2 


X2 + -> fxhXl -> ^"-"fc /i_/2 
X~2~ -» /l/ 2 Xl° -» /l / 2 


152 



TABLE XV: Number of events from gaugino pair production for Case(l) input at LHC (Tevatron benchmark (1)) with the 
final state involving t~t~&& + X (and charge conjugate channels). 



Channel 
number 


Gaugino 
pair 


a (pb) at LHC 
(Tevatron) 


Decay channel of 
the gauginos 


No of events 
at LHC (Tevatron) 


1. 


X2 U X2 U 


2.1 x 10~ 2 


X2° - / /Xi" -> I'cb f f 
X2° - / /Xi ° -> liCb / / 


54 


2. 


X2°Xl + 


5.3 x 10" 1 
(2.9 xlO^ 1 ) 


X2° - / /Xi" - 17 <*> f f 
Xi + -> /1 /2Xi° - I'cb f 1 h 


1500 
(73) 


3. 


Xi + Xi" 


4.7 x 10 _1 
(3.2 xlO^ 1 ) 


Xi + -» /1 /2Xi" ^ITcbf.h 
Xi -» /1 / 2 Xi° -> 'i c6 /1 / 2 


1318 

(85) 


4. 


X2 + X2~ 


2.3 x 10~ 2 


X2 + -/i hxi'^Kcb_f 1 / 2 

X2 ~> /l / 2 Xl° -> 'i C& /l / 2 


64 



TABLE XVI: Same as Table XV, but for Case (2) input, and benchmark (1) for Tevatron. 



VII. CONCLUSION 

In conclusion we reiterate that three relatively large 
input couplings of A^ fe type in the flavor basis at Mq 
each bearing a different lepton index can generate at the 
weak scale via RG evolution several additional A'-type 
couplings in the same basis (see Section III) with natu- 
rally suppressed magnitudes. This occurs due to the fla- 
vor violation inevitable in any model due to quark mix- 
ing. Also generated are three lepton number violating 
bilinear terms Ki . These bilinears and three induced A - 33 
couplings can explain the v oscillation data [18]. The 
upperbounds \' i33 < 1.5 x 10~ 4 and Ki < 1.2 x 1CT 4 
obtained from neutrino data can be translated into up- 
per bounds on the input couplings at Mq- In spite of 
the severe constraints from the v sector these input cou- 
plings are allowed to be reasonably large. When evolved 
to the weak scale and then rotated to the mass basis the 
corresponding physical couplings can trigger several in- 
teresting signals, which are not expected in models where 
the couplings in the v sector are the only ones with ap- 
preciable magnitude. 

In order to illustrate these ideas we study the RG equa- 
tions of the RPV mSUGRA model along with some mod- 
els of the quark mixing matrix [20, 25]. In models with 
quark mixing occuring in the up sector (Scenario I) we 



find that a minimal set of any three couplings, chosen 
from the group X' il3 and X' i23 , each with a different lepton 
index can serve the purpose. For the input set A- 13 (Case 
1), the oscillation data as well as other low energy con- 
straints [4, 16] lead to the upper bounds A' 113 < 0.0064, 
A 213 < 0.019, A 313 < 0.039 at M G for m ia = 100 GeV 
(see Sections III and IV). Rare r decays with BRs close 
to the current experimental bounds are allowed in this 
case (see Table I). For the input set A- 23 (Case 2), the 
bounds on the input couplings are more stringent, and 
hence no interesting rare weak process can be mediated 
by such small couplings. For the mixed choice A' 123 , A 213 
and A 313 , the rare decay K — > txvv with BR close to the 
current experimental upper limit are consistent with the 
bounds on the input couplings. 

LSP decays with characteristic BRs occur for all 
choices of GUT scale inputs. For Case 1 LSP decays into 
a charged lepton accompanied by a 6-jet and a light quark 
jet is the distinctive signal (Table II) . In Case 2 two 
heavy flavor jets (c and b) accompany the charged lep- 
ton (Table III) . The relative abundance of the final states 
with different leptons can indicate the relative strengths 
of the inputs at Mq. 

The lighter top squark decays in RPV channels can 
naturally compete with its RPC decays in models of v 
mass provided this sparticle happens to be the NLSP 
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[8, 9]. In Case 1, due to the relatively large upper bounds 
on the input RPV couplings, RPV can overwhelm RPC 
decays in a large region of the parameter space. However, 
if the input couplings (in particular the one involving r) 
are somewhat smaller than the upper bounds, competi- 
tion between RPV and RPC modes may show up (Table 
VII). Competition among the RPV and RPC decays are 
more probable in Case 2 (Table VIII). 

In Case 1 the upper bounds on the RPV couplings 
are consistent with direct RPV decays of other squarks 
and leptons with sizable BRs (Tables IV- VI). The above 
decays can lead to several low background signals both 
at the Tevatron and at the LHC. At Tevatron gaugino 
pair production followed by the LSP decay can lead to 
interesting signals (Tables XV and XVI). 

Top squark pair production at Tevatron followed by 
direct RPV decays of both the sparticles into an OSD 
pair of same flavor and jets have been discussed in [8, 9]. 
In this paper we consider the signals with isolated LSD 
pairs and isolated dileptons of different flavors (Tables IX 
and X corresponding to Case 1 and Case 2, respectively). 
Both are expected to be low background signals and our 
numerical estimates indicate large number of events. Ob- 
viously a huge signal at LHC is anticipated. 

Gaugino pair production followed by LSP decays can 
provide the signals similar to the ones discussed in the 
last paragraph. Even at the Tevatron a healthy size of 
the signal is expected in both Case 1 and 2 (Tables XIII 
and XIV). The prospect of observing these signals at the 
LHC is even better (Tables XV and XVI). 

Interesting signals arising from sfermion-antisfermion 
pair production followed by the direct RPV decay of one 
and indirect RPV decay of the other can lead to signals 



with decent magnitudes at LHC (Tables XI and XII). 

If, on the other hand, a model of the CKM matrix 
with mixing in the down quark sector (Scenario II) is 
considered, an mSUGRA type model with phennomcno- 
logically interesting, large input RPV couplings at Mq 
not directly related to the v sector is not viable. If at- 
tention is paid to the v constraint the input couplings at 
Mq are destined to be highly suppressed (Section VI) 

We next consider a phenomenological model with rel- 
atively large weak scale RPV parameters in the flavor 
basis induced by some high scale physics, the nature of 
which is presently unknown. Even if the induced trilinear 
couplings are not directly related to the neutrino sector, 
the A -33 couplings can be generated in the mass basis via 
CKM rotations. One can start with either (i) A^ 13 , i31 N 
or (ii) A- 23 ( i32 ) to be nonzero in the flavor basis at the 
weak scale. Unfortunately CKM rotations of these in- 
puts not only generate the couplings required by the v 
sector consistent with their upper limits, but also gener- 
ate other dangerously large couplings leading to unacc- 
etably large K° — K° mixing amplitude. Consequently 
such inputs must have small magnitudes not interesting 
for rare weak decays or direct RPV decays of squarks or 
sleptons. However top squark and LSP decays occur with 
characteristics quite different from the ones discussed for 
Scenario I (Section VI). 



VIII. ACKNOWLEDGEMENTS 

AK was supported by the project no. SR/S2/HEP- 
15/2003 of DST, Govt, of India. 



[1] Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998) (at- 
mospheric), S. Fukuda et al, Phys. Rev. Lett. 86, 5651 

(2001) , Q.R. Ahmed et al, Phys. Rev. Lett. 89, 011301 

(2002) , Phys. Rev. Lett. 89, 011302 (2002) (solar), M. 
Apollonio et al, Eur. Phys. J. C 27, 331 (2003), K. 
Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003) (re- 
actor), M.H. Ahn, Phys. Rev. Lett. 90, 041801 (2003) 
(long-baseline). 

[2] M. Gell-Mann, P. Ramond, and R. Slansky, in Supergrav- 
ity, ed. D. Freedman and P. van Nieuwenhuizen (North- 
Holland, Amsterdam, 1979), p. 315; T. Yanagida, in 
Proc. of the Workshop on Unified Theory and Baryon 
Number in the Universe, ed. O. Sawada and A. Sugamoto 
(KEK, Japan, 1979); R. Mohapatra and G. Senjanovic, 
Phys. Rev. Lett. 44, 912 (1980), Phys. Rev. D 23, 165 
(1981). 

[3] Some excellent reviews and text books are: H.P. Nilles, 
Phys. Rep. 110, 1 (1984); H.E. Haber and G.L. Kane, 
Phys. Rep. 117, 75 (1985); J. Wess and J. Bagger, Super- 
symmetry and Supergravity, 2nd ed., (Princeton, 1991); 
K. Intriligator and N. Seiberg, hep-ph/9509066; J.D. 
Lykken, hep-ph/9612114; S.P. Martin, hep-ph/9709356; 
N. Polonsky, hep-ph/0108236. M. Drees, R.M. Godbole, 



and P. Roy, Theory and Phenomenology of Sparticles 
(World Scientific, 2005). 

[4] For reviews on R-parity violating SUSY couplings and 
original references, see, e.g., H.K. Dreiner, in Perspectives 
on Super symmetry, ed. G.L. Kane, World Scientific (hep- 
ph/9707435); M. Chemtob, hep-ph/0406029; R. Barbier 
et al, hep-ph/0406039. 

[5] C.S. Aulakh and R.N. Mohapatra, Phys. Lett. B119, 136 
(1982); L.J. Hall and M. Suzuki, Nucl. Phys. B231, 419 
(1984); J.R. Ellis et al, Phys. Lett. B150, 142 (1985); 
GG. Ross and J.W. Valle, Phys. Lett. B151, 375 (1985); 
S. Dawson, Nucl. Phys. B261, 297 (1985); R. Barbieri et 
al, Phys. Lett. B252, 251 (1990); A.S. Joshipura and M. 
Nowakowski, Phys. Rev. D 51, 2421 (1995); F.M. Borzu- 
mati et al, Phys. Lett. B384, 123 (1996); Y. Grossman 
and H.E. Haber, Phys. Rev. Lett. 78, 3438 (1997), Phys. 
Rev. D 59, 093008 (1999); S. Davidson and M. Losada, J. 
High Energy Physics 0005, 021 (2000). Also see S. Rak- 
shit, Mod. Phys. Lett. A19, 2239 (2004), and reference 
[15] therein. 

[6] A. Datta and B. Mukhopadhyaya, Phys. Rev. Lett. 85, 
248 (2000); D. Restrepo, W. Porod and J.W.F. Valle, 
Phys. Rev. D 64, 055011 (2001). 



19 



[7] 
[8] 
[9] 
[10] 



[11] 
[12] 



[13] 

[14] 
[15] 
[16] 



[17] 

[18] 
[19] 



D. Acosta et al. (CDF collaboration), Phys. Rev. Lett. 
92, 051803 (2004). 

S. Chakrabarti, M. Guchait and N.K. Mondal, Phys. Rev. 
D 68, 015005 (2003); Phys. Lett. B600, 231 (2004). 
S.P. Das, A. Datta and M. Guchait, Phys. Rev. D 70, 
015009 (2004). [20 
R. Barbieri and A. Masiero, Nucl. Phys. B267, 679 ( 
1986); K. Agashe and M. Graesser, Phys. Rev. D 54, 4445 [21 

(1996) ; S.A. Abel, Phys. Lett. B410, 173 (1997); R. Bar- 
bieri, A. Strumia and Z. Berezhiani, Phys. Lett. B407, 

250 (1997); J.-H. Jang, Y.G. Kim and J.S. Lee, Phys. [22 
Lett. B408, 367 (1997); Phys. Rev. D 58, 035006 (1998); [23 
J.-H. Jang, J.K. Kim and J.S. Lee, Phys. Rev. D 55, 
7296 (1997); K. Huitu et al, Phys. Rev. Lett. 81, 4313 [24 

(1998) ; D. Choudhury, B. Dutta and A. Kundu, Phys. 
Lett. B456, 185 (1999); D. Chakraverty and D. Choud- 
hury, Phys. Rev. D 63, 075009 (2001); D. Chakraverty 
and D. Choudhury, Phys. Rev. D 63, 112002 (2001); H. [25 
Dreiner, G. Polesello and M. Thormeier, Phys. Rev. D [26 
65, 115006 (2002); A. Akeroyd and S. Recksiegel, Phys. 
Lett. B541, 121 (2002). [27 
J. P. Saha and A. Kundu, Phys. Rev. D 66, 054021 (2002). [28 

D. Guetta, Phys. Rev. D 58, 116008 (1998); G. Bhat- 
tacharyya and A. Datta, Phys. Rev. Lett. 83, 2300 [29 

(1999) ; G. Bhattacharyya, A. Datta and A. Kundu, Phys. 
Lett. B514, 47 (2001); A. Datta, Phys. Rev. D 66, 
071702 (2002); B. Dutta, C.S. Kim and S. Oh, Phys. [30 
Rev. Lett. 90, 011801 (2003); A. Kundu and T. Mi- 

tra, Phys. Rev. D 67, 116005 (2003); B. Dutta et al, 
hep-ph/0312388, hep-ph/0312389; G. Bhattacharyya, A. [31 
Datta, and A. Kundu, J. Phys. G 30, 1947 (2004). [32 
G. Bhattacharyya and A. Raychaudhuri, Phys. Rev. D [33 
57, 3837 (1998); [34 
J. P. Saha and A. Kundu, Phys. Rev. D 69, 016004 (2004). 
A. Kundu and J.P. Saha, Phys. Rev. D 70, 096002 (2004). 
V. Barger, G.F. Giudice and T. Han, Phys. Rev. D 40, 
2987 (1989); G. Bhattacharyya, hep-ph/9709395; B.C. [35 
Allanach, A. Dedes and H.K. Dreiner, Phys. Rev. D D60, 
075014 (1999). Also see M. Chemtob, and R. Barbier et [36 
al, in [4]. 

M. Tegmark et al. (SDSS Collaboration), Phys. Rev. D 
69, 103501 (2004). [37 
A. Abada and M. Losada, Phys. Lett. B492, 310 (2000). 
S. Roy and B. Mukhopadhyaya, Phys. Rev. D 55, 7020 [38 

(1997) ; M. Drees et al, Phys. Rev. D 57, 5335 (1998); 

E. J. Chun et al, Nucl. Phys. B 544, 89 (1999); S. 
Rakshit, G. Bhattacharyya and A. Raychaudhuri, Phys. 



Rev. D 59, 091701 (1999); M. Hirsch et al, Phys. Lett. 
B486, 255 (2000); A. Datta, B. Mukhopadhyaya and S. 
Roy, Phys. Rev. D 61, 055006 (2000); A.S. Joshipura, 
R.D. Vaidya, and S.K. Vempati, Nucl. Phys. B639, 290 
(2002). 

B.C. Allanach, A. Dedes and H.K. Dreiner, Phys. Rev. 
D D69, 115002 (2004). 

R. Hempfling, Nucl. Phys. B478, 3 (1996); M. Hirsch et 
al, Phys. Rev. D 62, 113008 (2000), Phys. Rev. D 65, 
119901(E) (2002). 

F. de Campos et al, Phys. Rev. D 71, 075001 (2005). 

B. de Carlos and P.L. White, Phys. Rev. D 54, 3427 

(1996), Phys. Rev. D 55, 4222 (1997). 

L.J. Hall and M. Suzuki in [5]; D.E. Brahm and L.J. Hall, 

Phys. Rev. D 40, 2449 (1989); K.Tamvakis, Phys. Lett. 

B382, 251 (1996); G.F. Giudice and R. Rattazzi, Phys. 

Lett. B406, 321 (1997). 

K. Agashe and M. Graesser, in [10]. 

S. Eidelman et al. (Particle Data Group Collaboration), 

Phys. Lett. B592, 1 (2004). 

F.E. Paige et al, hep-ph/0312045. 

L.E. Ibanez and G.G. Ross, Phys. Lett. B110, 215 
(1982). 

B. Gato, Nucl. Phys. B 278, 189 (1986); P. Moxhay and 
K. Yamamoto, Nucl. Phys. B 256, 130 (1995); N. Polon- 
sky and A. Pomarol, Phys. Rev. D 51, 6532 (1995). 

M. Drees, Phys. Lett. B181, 279 (1986); Y. Kawamura, 
H. Murayama and M. Yamaguchi, Phys. Rev. D 51, 1337 
(1995). 

S.P. Das, A. Datta, and S. Poddar, work in progress. 
V. Barger et al, Phys. Lett. B538, 346 (2002). 
http: //theory . sinp.msu.ru/~pukhov/calchep .html 
See, e.g., R. Yamada, Fermilab.FN.629, D0note 2250 
(1994); R.J. Lipton for the CDF and D0 collaborations, 
hep-ex/0505075; S. D'Auria for the CDF and D0 collab- 
orations, hep-ex/0506049. 

M. Drees, Phys. Lett. B158, 409 (1985); J. Ellis et al, 
Phys. Lett. B155, 381 (1985). 

A. Datta et al, Phys. Rev. D 61, 055003 (2000); A. 
Datta, A. Datta, and S. Maity, J. Phys. G 27, 1547 
(2001). 

K.I. Hikasa and M. Kobayashi, Phys. Rev. D 36, 724 
(1987). 

C. Boehm, A. Djouadi, and Y. Mambrini, Phys. Rev. D 
61, 095006 (2000). 



